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1. Introduction 

Materials with low optical absorption coefficients are essential for high-precision laser- 
interferometric measurements. Absorptions in mirror substrates of as low as 10~^/cm al- 
ready limit gravitational wave detectors because absorption leads to heating and a thermal 
deformation of the mirrors [1] and also to photo-thermal noise [2]. Future gravitational wave 
detectors will use cryogenically cooled mirrors [3] to reduce thermally excited motions of 
mirror surfaces. Then, low optical absorptions will become even more crucial. Consequently, 



the measurement of small absorption coefficients in the regime below 10~^/cm is important 
to find appropriate mirror materials and to enable the reliable design of future gravitational 
wave detectors, such as the Einstein Telescope [4,5]. 

In the past, several methods have been developed that are able to measure absorption 
coefficients of the order of 10~^/cm. All these methods are based on indirect measurement 
schemes. They do not directly sense the power loss of a transmitted beam but utilize the 
temperature increase that arises due to the absorption. In calorimetric approaches the tem- 
perature increase is directly measured [6]. Other approaches exploit light beam defiection or 
beam shape deformation due to local heating [7, 8] . 

In this paper we present another indirect measurement scheme to determine small ab- 
sorptions. The material under investigation is put inside an optical cavity whose length is 
linearly scanned over a cavity Airy peak. Approaching cavity resonance the temperature 
along the cavity mode increases and the optical path length for a cavity round trip changes. 
The thermally induced optical path length change is a photo-thermal self-phase modulation 
resulting in a deformed shape of the Airy peak. Since the phase change depends on the light 
intensity it may be considered as the result of a "thermo-optic Kerr-effect" . Importantly, the 
Airy peak deformation depends on the scan direction, i.e. whether the cavity is shortened or 
lengthened. The hysteresis in the time-resolved measurements provides information of the 
absorption coefficient, if relevant material parameters are known and included in a numerical 
time-domain simulation. A positive side-effect of our method is the power build-up inside the 
cavity which compensates the need for laser sources with higher powers when approaching 
the regime of extremely low absorption. 

2. Theory and Method 

In this section we describe the time-domain simulation that is used to analyse the measure- 
ment data and to deduce the absorption from it. Our approach is based on work by Hello and 
Vinet [9, 10] in which they describe the heating of an absorbing material due to a Gaussian 
laser beam. In our case, a sample of the absorbing material with polished (plane) surfaces 
is placed inside a (high) finesse cavity. One may choose the light's angle of incidence to 
be the Brewster angle to avoid refiection losses. A schematic is shown in Fig. 1. When the 
cavity round-trip phase 0cav is linearly increased (or reduced) by 6{t) and scanned over a 
cavity resonance, absorption leads to a dynamic temperature profile inside the material and 
inside the cavity mirror surfaces. The result is a (photo-thermal) self-phase modulation and 
a deformation of the cavity Airy peak. 

Our time-domain model iteratively calculates the intracavity intensity after each round 
trip. The time t is discretized becoming an integer multiple of the round trip time, yielding 
t = tn = n/A/psR, where A/fsr is the cavity free-spectral-range. The intra-cavity field 



ttn := a{tn) after n round trips reads 



a„ = i^l - rl aine'-^'-^*") + nh e'"^"^'"^ an-i • (1) 

Here, ri is the amplitude reflectivity of the first mirror, whereas r2 is the effective amplitude 
reflectivity of the second mirror, which includes all round trip losses. The amplitude of the 
incident power Pjq is given by a-m = 2^^/I^^JJeocyi^w^, where wq is the waist radius of the 
beam eo the dielectric constant and c the speed of light. The cavity input field gains the 
phase (pin that is due to the temperature gradient inside the incoupling mirror emerging 
from its coating absorption. The phase 0„(t„) after n round trips can be written as 

4>n = S(tn) + 0spm(^n., «) , 

where (5(t„,) is the phase due to the external cavity detuning and 0spm(^n5«) is due to the 
photo-thermal (internal) self-phase-modulation which depends on the absorption a. The 
external detuning for the round trip number n is determined from 

S{tn) = 60 + n27rNFSR-r^ ■ (2) 

^JFSR 

Here, A^fsr is the number of free spectral ranges that were scanned with frequency Ug- The 
velocity v^i of the scanning mirror is therefore given by 

Vm = ^XUs ■ NpsR ■ (3) 

The temperature distribution T„ for round trip n is calculated by using the recurrence 
relations (Eq. (15) in [10]). These equations determine the radial and longitudinal tempera- 
ture gradient at any time including thermal conductivity. The starting point is the external 
temperature Tq. The detuning 0spm induced by the photo-thermal self-phase- modulation is 
then given by equations (33) and (35) of [9]. Note that 0spm includes the effects of a non-zero 
thermo-optic coefficient dn/dT and a non-zero expansion coefficient Oth- As a starting point 
for the numerical simulation we use the steady-state solution for the start detuning Sq 



Assuming a perfect mode-matching of the input field Om to the cavity mode, the reflected 
and transmitted flelds are given by 



arcfi(^n) = iy 1 - rf a„(a) + naine"^'" , 

^ 

atrans(in) = iv 1 " r2^a.„(a) . 

Eq. (5) can now be used to calculate the time-resolved shape of an Airy peak. By varying 
a the result can be fltted to the measurement performed in reflection or transmission of 



the cavity. Figure 1(b) shows an example, i.e. simulated Airy peaks obtained from a cavity 
containing some absorbing material with dn/dT > 0. The resonance peaks get broader 
when shortening the cavity because the positive thermo-refractive coefficient counteracts the 
external change of the cavity length. Accordingly the resonance peaks get narrower when 
the cavity length is increased. In particular the hysteresis can be used to precisely determine 
the absorption of the material. For comparison we also show the normal Airy peak without 
self-phase modulation (dashed line). 

3. Measurements and data analysis 

To characterize the feasibility of our method we performed a series of absorption measure- 
ments on a 7mol% MgO-doped LiNbOa crystal. A single measurement set involves a charac- 
terization of the piezo electric element that is used to change the cavity length, and altogether 
four time-resolved photo-electric detections. A fast photo-diode records the Airy peaks in 
reflection of the cavity when the latter is (a) lengthened or (b) shortened, using (1) a low 
laser power without any thermal Airy peak deformation or (2) a laser power at which a 
thermal deformation is clearly visible. The low-power setting is used to quantify the two 
reflectivities rf and f|. The high-power setting is used to quantify the absorption coefficient 
ttLN- All three quantities and their error bars are deduced from a single measurement set 
and a numerical time-domain simulation applying a Markov-chain Monte-Carlo (MCMC) 
algorithm. Records (al) and (bl) are identical thereby con&ming that the laser power was 
chosen to be low enough so that thermal effects do not yet come into play. Note that (al) 
and (bl) will not necessarily have the shape of the central Airy peak (dashed) shown in 
Fig. lb but may show a ringing effect due to the cavity loading or decay time [21]. This 
effect is also precisely modelled in our simulation. 

3. A. Experimental setup 

In our research group, we routinely use cavities containing MgO-doped LiNbOa (LN) crystals 
for second harmonic generation (SHG) and squeezed light generation (SLG) at a wavelength 
of 1064 nm [19,22,23]. The optical absorption of these nonlinear crystals is a limiting factor 
in achieving high conversion efficiencies and high squeezing factors. Accurate absorption 
coefficients are therefore required to optimize the nonlinear cavity design. Unfortunately, 
manufacturers' data typically are rather inaccurate and a standard value of oln ^ 10^'^/cm 
at 1064 nm is quoted in most cases. In this work we used one of our SHG cavities to measure 
the absorption coefficient of LN and test our new absorption measurement technique. 

Figure 2 shows the experimental setup. An incoupling mirror and the curved side of the 
plano-convex LiNbOa crystal form a single-ended standing wave cavity for laser light at 
1064 nm. The cavity mirrors have power reflectivities of r^ = i?i ?a 90 % and rg = -R2 > 



99.8 %, respectively. A small air gap separates the incoupling mirror from the anti-reflection 
coated, plane crystal surface. Table 1 contains detailed geometric parameters of this resonator 
and the laser beam as well as the material parameters of the LiNbOa-crystal. 

Up to 1.5 W of single mode radiation at 1064 nm was modematched into the cavity with a 
modematching efficiency of greater than 95%. To prevent the generation of second harmonic 
radiation, both the input field polarization and the crystal temperature were detuned from 
their usual operation point. A piezoelectric transducer (PZT) moved the incoupling mirror 
to allow for a scan of the cavity length. The photo diode measured the temporal behaviour 
of the refiected laser power. We ensured that the photodiode was fast enough, i.e. had a high 
bandwidth, so that it did not influence the shape of the recorded Airy peaks. 

Fig. 3 shows an example of Airy peaks with visible thermal effects as measured in reflection 
of the cavity. The blue curve forms for a lengthening resonator, the red curve for a shortening 
resonator. No parameter other than the scan direction was changed. The two curves would 
be identical without self-phase modulation and no hysteresis effect would occur without 
absorption. The solid lines in Fig. 3 represent our simulation fltted to the experimental data. 
The narrow curves show a discrepancy in the left wings, the broad curves show a discrepancy 
in the right wings. The two deviations come from the non-perfect modematching to the cavity 
and the excitation of a higher-order cavity mode. 

Apart from taking a simple full measurement set, i.e. lengthened and shortened resonator 
at two different laser powers, we performed measurements at three different laser powers and 
three different scan frequencies. While not strictly necessary for an absorption measurement, 
these measurements demonstrate the consistency of our result, see below. 

3.B. Measurement analysis 

For the analysis of the measured peaks the PZT had to be calibrated because of its own 
hysteresis and non-linearity. This calibration was done at low laser powers where no thermal 
effect occurred. We measured the width of the Airy peaks at different positions of the PZT's 
scanning range by slightly shifting the laser frequency. A third-degree polynomial well de- 
scribed the peak width depending on peak position. Together with a scan showing a full free 
spectral range, we used this polynomial to linearize the PZT movement. We performed the 
calibration for both scan directions and for each scan velocity that we used. Measurements 
were performed at three different scan velocities, namely 2 ■ 1064 nm/5 ms, 2 ■ 1064 nm/2.5 ms 
and 2 ■ 1064 nm/0. 285 ms. For each scan velocity we measured Airy peaks at three different 
input powers: 100 mW, 750 mW and 1.5 W. 

Table 1 gives a complete list of the parameters that enter our simulation. We used values 
from literature for the material parameters, for the geometric parameters we chose values to 
our best knowledge of the cavity design. The mirror reflectivities Ri and i?2 deflne the cavity 



resonance width and power build-up. Here, R2 is an effective reffectivity, which includes 
absorption and scattering losses. Only this value, rather than the pure reffectivity R2, is 
accessible when light enters the cavity through mirror Ri. As resonance width and power 
build-up have a strong impact on the heating of the substrate, we do not use the reffectivity 
values as given by the coating manufacturer. Instead, we treat Ri and R2, as well as the 
absorption oln; as free parameters of our simulation. 

For low input powers and fast scan velocities the resulting temperature change inside the 
substrate is small and no deformation of the peaks is visible. Such time series are optimally 
suited to determine Ri and i?2- Towards higher input powers and lower scan velocities, 
the peaks begin to show a hysteresis. For all our measurements that were performed with 
different laser powers, the hysteresis values could be explained completely by the self-phase 
modulation, i.e. our simulation provided a very good discription of the measurement. From 
this we conclude that no spatial deformation of the cavity mode occurred. 

We performed a quantitative analysis by calculating the variance between simulated and 
measured data. Starting from an initial set of parameters, we ran a Metropolis-Hastings 
Markov-chain Monte Carlo (MCMC) [20] algorithm which minimized the variance. The data 
chains generated can be converted into histograms for the free simulation parameters. The 
histograms for the reflectivities -Ri, R2 and for the absorption ojln as derived from a single 
measurement setting are shown in Figure 4. As the histograms closely resemble Gaussian 
distributions, we give the mean value and standard deviation of all nine measurements in 
Table 2. 

The mean value of the results for the incoupling mirror was found to be -Ri = (89.43 ± 
0.75) % which is in good agreement with the manufacturer's specification for this coating 
((90 ± 1) %). The effective value for the high-reflective coating of the crystal was determined 
to be R2 = (99.79±0.01) %, which is also in accordance with the specifications. The measure- 
ment with 100 mW input-power at a scan velocity of f = 2 • 1064 nm/5 ms was the boundary 
where a small thermal effect was visible. However, no accurate absorption coefficient could 
be deduced due to rather large error bars. Four measurements showed a significant thermal 
effect and were used to derive four independent values for the absorption coefficient of LN. 
All four values for oln have mutually overlapping error bars. Figure 5 gives a graphical 
overview of the results for oln for all measurements. The mean value of the four results is 
«LN = 5.9 X 10~^/cm. As the error bar we quote the standard deviation (an averaged value) 
of the single measurement set which typically was ±0.8 x 10~^/cm. This number includes 
the influence from errors in the reflectivities Ri and R2 as shown in Fig. 4. 



3. C. Error propagation 

We considered the influence of possible errors in tlie input parameters on the resulting value 
for oln (from a single measurement). For this investigation we individually changed the 
values of the simulation input parameters and recalculated i?i, i?2 and oln for each case. 

Our investigation showed that the parameters can be grouped into two categories. The first 
category contains parameters that have a very weak influence on the absorption coefficient in 
our case. For our system, heat radiation described by the material emissivity 0.0 < e < 1.0 
is not relevant at all, because the substrate is heated only within the beam radius, far 
away from the substrate's surface. Also the absorption coefficient of the substrate coatings 
(ttcoating) can be neglected, as it is much smaller than the substrate absorption oln and the 
coating thickness is negligible compared to the substrate dimension. A few percent change 
of the values for the index of refraction n, the intra-cavity airgap s, the substrate radius 
R, and the beam waist uq also has a negligible effect on the absorption coefficient oln- 
The second category contains the remaining parameters of our model. These parameters and 
their respective influence on a^N for a 4 % change in the parameter value are the input laser 
power P (3.6%), the substrate length L (3.6%), the thermal conductivity /cth (1.6%), the 
thermal refractive coefficient dn/dT (3.6%), the thermal expansion Oth (1-4%), the density 
p and the heat capacity c (2.8%). Note that in the simulation p and c always appear as a 
product, and the influence of their error bars is identical. 

Assuming that our measured parameters as well as the material parameters from literature 
are precise to within 4 % and statistically independent from each other, we conclude that the 
error of ±0.8 x 10~^/cm (±13.6 %) coming directly out of the Markov-chain Monte-Carlo 
simulation dominates the error on our final result. The total error sums up to ±0.9 x 10~^/cm 
(±15.7%). 

3.D. Sensitivity of the method 

To make a prediction of the sensitivity of our method, we consider the absorption measure- 
ment of crystalline silicon at a wavelength of 1550 nm. This value has not been measured 
before, but data at shorter wavelengths [24,25] suggest an absorption coefficient smaller than 
10~^/cm in case of pure silicon [26]. Our simulation is based on a 6.5 cm long silicon sample 
inside a cavity of finesse 20,000 pumped with 1 W of input laser power. The reflectivities of 
the two cavity mirrors are assumed to be identical. Fig. 6 a shows the Airy peaks as detected 
in the reflected light for both cavity scan directions with an absorption of 10~^/cm. The 
scan velocity used in that simulation was v^ = 2 ■ 1550 nm/s. The curves are normalized to 
the input power of 1 W. Both curves show oscillations and values above unity which arise 
from the cavity loading and decay time [21]. Fig. 6 b shows the difference of the two scan 
directions normalized to the Airy peak without absorption. We find a significant hysteresis 
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curve that reaches up to 12% of the input power. Our simulation neglects the influence of 
the absorption in the dielectric coatings. In practice, the absorption inside the cavity mirror 
coatings has to be insignificant as in our experiment or it has to be measured independently 
when the sample is removed from the cavity. Anti-reflection coatings or high-reflection coat- 
ings on the sample itself can also be taken into account when two different sample lengths 
are studied. Generally, the photo-thermal self-phase modulation from the coating absorption 
must not dominate the overall photo-thermal effect inside the cavity. New low-loss coating 
materials such as diamond [27] or monolithic, nano-structured surfaces [28] might be used. 

4. Conclusion 

In this paper we introduce a new low-absorption measurement method based on the optical 
phase change inside the material when absorption leads to local heating. The effect is 
understood as a cavity- assisted photo-thermal self-phase modulation of light. We used our 
method to determine the absorption coefficient oln of a LiNbOs crystal. Our result of 
QiLN = (5.9 ± 0.9) X 10~'^/cm is in accordance with the typically referred upper bound of 
10~^/cm as available on manufacturer websites. Measurements with different laser powers 
could all be well described without considering a spatial mode distortion. We conclude that 
no such mode conversion occurred in our experiments. However, this might be possible 
at even higher laser powers or smaller waist sizes. We theoretically applied our method 
to a material with an absorption coefficient of a = 10~^/cm. We conclude that such low 
absorptions should be measureable when a sample of a few cm length is put into a cavity 
with a finesse of the order 10^. Our time-resolved Markov-chain Monte-Carlo (MCMC) 
simulation is based on a variety of material parameters. The coupling of parameter errors 
into the error of the absorption coefficient is linear or less. 
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List of Figures 

1 Scheme of the measurement method, (a) When changing the cavity round trip 
phase by S{t), absorption of intracavity power leads to a temperature change 
T{t) and an additional phase for the intracavity field (pcav Ri and R2 are 
the mirror power refiectivities; 0in is the input phase. By using the Brewster 
angle, surface refiections can be avoided, (b) Airy peaks showing hysteresis 
due to the photo-thermal self-phase modulation. The dashed line shows the 
peak without any absorption {a = 0), whereas the blue (narrow) and red 
(broad) peaks show the Airy peak for the same absorption coefficient a > 0, 

but for lengthening and shortening the cavity, respectively 14 

2 Experimental setup: The resonator is formed by the incoupling mirror with 
reflectivity -Ri ~ 90 % and the crystal's highly reflecting (HR) coating with a 
reflectivity of R2 > 99.8 %. The resonator length is scanned with a frequency 
from a function generator (FG) which is fed through a high voltage amplifier 
(HV) 15 

3 Example of measured (solid) and simulated (dashed) Airy peaks with visible 
thermal effect. Without absorption all curves would be identical. The curves 
were measured in reffection, no parameter other than the scan direction was 
changed. The red curve forms for a shortening resonator, the blue one for a 
lengthening resonator 15 

4 The Metropolis-Hastings MCMC algorithm draws samples from the simula- 
tion parameter space creating a chain of individual realizations that result in 
the parameter distributions. Here we show the histograms of a chain obtained 
from a single measurement set at a laser input power of 0.75 W and a scan- 
velocity oi V = 2 ■ 1550 nm/2. 5 ms. -Ri (top) and R2 (middle) are required to 
characterize the cavity. The bottom figure shows the result for oln- The bars 
represent histograms of the MCMC run. The curves are gaussian fits to the 
histogramms 16 

5 Four independent measurement values of the absorption coefficient oln and 
their statistical standard deviations. The blue lines show the mean value of 
the four measurements, which is oln = 5.9 x 10~^/cm. The dashed blue lines 
mark the averaged standard deviation of Aq^ln = 0.8 x 10~^/cm 17 
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Simulated hysteresis effect for the Airy peaks in reflection from a monolithic 
silicon cavity of finesse 20000. The curves are normalized to the input power 
of 1 W at 1550 nm. The scan velocity of the cavity length is 2 • 1550 nm/s and 
the absorption was assumed to be 10~^/cm. Other parameters can be found in 
Table 1. (a) Airy peaks for lengthening (blue) and shortening (red) the cavity, 
(b) The difference of the two scan directions AP normalized to the incident 
laser power of 1 W 1^ 
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Table 1: Material and geometric parameters of the LiNbOa- and Si-samples and cavity geo- 
metric parameters used for the simulations. 



Material parameters 


LiNbOs 




Si 


index of refraction n 


2.147 [11 




3.48 [15] 


thermal refr. coeff. dra/dT 


38.5 • 10-' 


VK [17] 


176.0 -lO-e/K [15 


specific heat c 


630J/(kgK) [13 


713J/(kgK) [14 


density p 


4635 kg/m^ [12; 


2330 kg/m^ [13 


thermal expansion Oth 


14.8 ■ 10-' 


VK [13] 


2.53-10-VK [14 


thermal conductivity k^h 


4.19W/(mK) [18 


1.56W/(mK) [16 


material emissivity e 


1.0" 




1.0" 


coating absorption acoatmg 


0.0 /cm 




0.0 /cm 


Cavity geometric parameters 


airgap s 


24 mm 




mm 


beam waist ujq 


24 /im 




160 //.m 


crystal length L 


6.5 mm 




65.0 mm 


crystal radius R 


2 mm 




50.0 mm 



"0.0 < e < 1.0 are the boundaries for the thermal emissivity. For our systems the value of 
this parameter is not relevant since R^ Uq. 
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Fig. 1: Scheme of the measurement method, (a) When changing the cavity round trip phase 
by S(t), absorption of intracavity power leads to a temperature change T(t) and an additional 
phase for the intracavity field 0cav Ri and R2 are the mirror power reflectivities; 0in is the 
input phase. By using the Brewster angle, surface reflections can be avoided, (b) Airy peaks 
showing hysteresis due to the photo-thermal self-phase modulation. The dashed line shows 
the peak without any absorption (a = 0), whereas the blue (narrow) and red (broad) peaks 
show the Airy peak for the same absorption coefficient a > 0, but for lengthening and 
shortening the cavity, respectively. 
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Fig. 2: Experimental setup: The resonator is formed by the incouphng mirror with reflectivity 
-Ri ~ 90 % and the crystal's highly reflecting (HR) coating with a reflectivity of i?2 > 99.8 %. 
The resonator length is scanned with a frequency from a function generator (FG) which is 
fed through a high voltage amplifier (HV). 
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Fig. 3: Example of measured (solid) and simulated (dashed) Airy peaks with visible ther- 
mal effect. Without absorption all curves would be identical. The curves were measured in 
reflection, no parameter other than the scan direction was changed. The red curve forms for 
a shortening resonator, the blue one for a lengthening resonator. 
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Fig. 4: The Metropolis-Hastings MCMC algorithm draws samples from the simulation pa- 
rameter space creating a chain of individual realizations that result in the parameter dis- 
tributions. Here we show the histograms of a chain obtained from a single measurement set 
at a laser input power of 0.75 W and a scan-velocity of t' = 2 ■ 1550 nm/2.5 ms. -Ri (top) and 
i?2 (middle) are required to characterize the cavity. The bottom figure shows the result for 
ctLN- The bars represent histograms of the MCMC run. The curves are gaussian fits to the 
histogramms. 
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Fig. 5: Four independent measurement values of the absorption coefficient oln and their 
statistical standard deviations. The blue lines show the mean value of the four measurements, 
which is oln = 5.9 x 10~^/cm. The dashed blue lines mark the averaged standard deviation 
of AttLN = 0.8 X lO^Vcm. 



Table 2: Results for i?i, R^ and oln: Mean values as well as standard deviation of the 
parameters are given. 
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Fig. 6: Simulated hysteresis effect for the Airy peaks in reflection from a monolithic silicon 
cavity of finesse 20000. The curves are normalized to the input power of 1 W at 1550 nm. 
The scan velocity of the cavity length is 2 • 1550 nm/s and the absorption was assumed to be 
10~^/cm. Other parameters can be found in Table 1. (a) Airy peaks for lengthening (blue) 
and shortening (red) the cavity, (b) The difference of the two scan directions AP normalized 
to the incident laser power of 1 W. 



